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INTRODUCTION 


The interest in advanced tactical aircraft designed for efficient cruise and 
maneuver at supersonic speeds has highlighted the limitations of supersonic linear 
theory. The linear theory is well suited to slender transport configurations which 
satisfy the thin-wing and small-disturbance assumptions of the method. However, at 
supersonic speeds, the tactical-aircraft characteristics of low fineness ratio, 
rounded-wing leading edges, and moderate wing sweep, which result in transonic Mach 
numbers normal to the wing leading edge, present a formidable challenge to the 
linear-theory methods. Perhaps the most demanding problem occurs for the high-lift 
conditions required for supersonic maneuver. 

Basically, two approaches are available for the design of wings to produce low 
drag due to lift at high-lift conditions. One approach, which has been demonstrated 
experimentally at subsonic speeds, is to use a sharp leading-edge flap to produce 
separated flow and maintain a leading-edge vortex which provides vortex lift and some 
effective leading-edge suction. The second approach is to provide an attached-flow, 
controlled expansion around the leading edge of the wing. This latter approach is 
the subject of this report. 

To produce high lift with attached flow at supersonic speeds, the flow must 
accelerate to conditions at which the cross-flow velocity is supercritical. The 
basic idea is to generate high levels of lift using the low pressures resulting from 
the upper-surface supercritical cross flow while minimizing drag by avoiding large 
pressure gradients which separate the flow and by avoiding strong shocks which result 
in energy losses. The concept of controlling this supercritical cross flow at 
supersonic speeds (ref. 1 ) is a natural extension of the well-understood concepts 
developed for supercritical airfoils at transonic speeds. 

In order to accurately analyze and/or design wings with supercritical cross 
flow, it was necessary to have a computer code capable of accurately and effi- 
ciently analyzing highly nonlinear supersonic flows. To meet these requirements, the 
development of a series of full-potential supersonic flow codes (refs. 2 to 6) has 
been an integral part of developing the wing-design concept. Initially, a conically 
cambered wing was designed using the conical nonlinear potential code. This conical- 
wing experiment proved that the high-lift, supercritical-cross-flow wing-design con- 
cept was valid and that the recompression of the supercritical cross flow could be 
controlled to avoid boundary-layer separation (refs. 7 and 8). Subsequently, a 
three-dimensional cambered wing representative of wing planforms resulting from 
advanced tactical-fighter studies (ref. 9) was designed using the three-dimensional 
nonlinear full-potential code (NCOREL, ref. 6). 

The purpose of this paper is to present results of the experimental validation 
for the three-dimensional cambered wing which was designed to achieve attached 
supercritical cross flow for lifting conditions typical of supersonic maneuver. The 
design point was a lift coefficient of 0.4 at Mach 1.62 and 12° angle of attack. 
Results from the nonlinear full-potential method are presented to show the validity 
of the design process along with results from linear-theory codes. Longitudinal 
force and moment data and static-pressure data were obtained in the Langley Unitary 
Plan Wind Tunnel (ref. 10) at Mach numbers of 1.58, 1.62, 1.66, 1.70, and 2.00 over 


an angle-of -attack range of 0° to 14° at a Reynolds number of 2.0 x 10^ per foot. 
Oil-flow photographs of the upper surface were obtained at M = 1.62 for a « 8°, 
10°, 12°, and 14°. 


SYMBOLS 

The moment reference point is 16.701 in. behind the model apex on the centerline 
and 0.275 in. below the model reference line. Symbols in parentheses are used in 
some appendix tables and figures. 
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root chord length, 23.84 in. 
spherical marching-step size in NCOREL 
lift-drag ratio 
leading edge 
free-stream Mach number 

cross -flow Mach number, 

Mach number normal to leading edge, M cos A 
local static pressure 
free-stream stagnation pressure 
free-stream static pressure 
free-stream dynamic pressure 
free-stream Reynolds number, per foot 
wing leading-edge radius 

o 

reference wing area, 342.11 in 

free-stream stagnation temperature 

lateral perturbation velocity component 

vertical perturbation velocity component 

longitudinal distance measured from model apex, in. 

spanwise distance measured from model centerline, in. 

vertical distance measured from model reference plane, 
positive up, in. 

angle of attack, deg 

angle of attack at zero lift, deg 

angle between horizontal and circular-arc camber line at 
wing leading edge (see fig. 4) 
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local nondimensionalxzed spanwise coordinate, — — 
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ORIGINAL PAGE \9 
OF POOR QUALITY 


0 T streamwise airfoil twist angle, deg (see fig. 5) 

A leading-edge sweep angle, deg 

Subscripts : 

scp supercritical panel 

LE leading edge 

TE trailing edge 


AERODYNAMIC DESIGN 

The left-hand side of the following sketch illustrates the typical high -lift 
pressure distribution on an uncambered spanwise wing section with a rounded leading 
edge. The right-hand side shows the desirable pressure distribution of a properly 
shaped wing section. The proper camber and thickness eliminates the leading-edge 
expansion peak and reduces the strength of the cross-flow shock. The resultant 
upper-surface pressure distribution features both a supercritical cross-flow region 
(M c > 1) and a subcritical cross-flow region (M < 1). Tine attached supercritical- 
cross-flow concept attempts to maintain attached flow so that the drag reduction 


LEADING-EDGE 
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created by the pressure expansion on the rounded leading edge and on the forward- 
facing upper-surface slopes can be used to improve the wing performance. 

Design conditions of M = 1.62 and = 0. 4 were chosen as representative of 
future tactical-aircraft maneuver conditions. The wing planform selected for inves- 
tigation (fig. 1) was derived from an advanced tactical-fighter study (ref. 9). The 
basic leading-edge sweep angle was 57°, which corresponds to b cot G = 0.83 and 
M n = 0.88 at the design Mach number. Inboard of about 44 percent semispan, the wing 
was blended into a 65° leading-edge sweep angle. The outboard trailing-edge sweep 
angle was 33°, which blended into an 11° trailing-edge sweep angle inboard of about 
54 percent semispan. 

Given the design conditions and wing planform, the aerodynamic design problem is 
to specify a target pressure distribution and define a wing camber and thickness 
shape which generates the target pressure distribution. To aid in obtaining a target 
pressure distribution, a procedure for assessing the effect of variations in the size 
of the supercritical cross-flow region and in the pressure level for that supercrit- 
ical cross-flow region was developed using a modified linear-theory code described in 
reference 11. This procedure allows the specification of a conical region of super- 
critical cross flow of arbitrary size and pressure level near the wing leading edge. 
In the presence of these supercritical panel pressures, the subcritical wing pres- 
sures are then determined to minimize the drag due to lift of the entire wing for the 
specified Mach number and lift coefficient. The results of a typical design exercise 
are shown in figure 2. Each curve in the figure represents a different size (denoted 
by hscp^ of t ^ ie su P ercr itical panel; the variation of drag due to lift is shown for 
supercritical-panel pressure levels ACp^g C p ranging from 0.42 to 0.52. The chosen 
target is less than 5 percent above the minimum drag level and represents a pressure 
level and size for the supercritical cross flow which is intuitively felt to be 
attainable in the real flow. 

The wing-design target pressure distribution must make the transition from 
supercritical to subcritical cross-flow conditions, and it is desirable for this 
transition to occur isentropically (shockless). If, however, shocks cannot be 
avoided, their strengths should be controlled to maintain' low wave drag and not 
separate the flow. According to two-dimensional experimental data summarized in 
reference 12, shocks which produce static-pressure increases of less than 50 percent 
will not cause flow separation? therefore, in this wing design, all shock -produced 
static-pressure jumps were kept below 25 percent. 

The specification of an airfoil to produce the target pressure distribution was 
accomplished by using the nonlinear flow method of reference 6 (NCOREL) to design by 
iteration. The computer code solves the supersonic, full-potential equation using 
the exact surface boundary conditions. Therefore, the method treats the surface 
shape instead of computing thickness and camber effects independently. It also pro- 
vides accurate information at the leading edge, which is in contrast to the state-of- 
the-art linear potential theory. As a means of simplifying the airfoil development, 
the thickness envelope is generated first, and then the camber surface is generated. 

The wing leading-edge geometry was found to be critically important in attached- 
flow, high-lift design. The leading-edge radius is required to be large, by conven- 
tional supersonic wing-design standards, to prevent flow separation on the highly 
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loaded leading edge. A modified NACA four-digit thickness distribution was selected 
because the leading-edge radius can be easily varied using analytic equations which 
define the thickness distribution. The airfoil thickness shape selected corresponds 
to a leading-edge radius distribution shown in figure 3 with the maximum thickness 
ratio of 4 percent located at 40 percent of local chord. 

Once the airfoil thickness envelope was established, a systematic means of 
developing the camber surface was employed. An analytic description of the wing was 
obtained by superimposing the following three basic camber elements: spanwise 

circular-arc camber, dihedral, and twist. These three basic camber elements were 
systematically varied to obtain the spanwise target pressure distribution at five 
longitudinal control stations; the control stations were arbitrarily selected to be 
at 5, 10, 15, 20, and 25 in. aft of the wing apex. This procedure resulted in a 
dihedral angle of 10°, a longitudinal distribution of circular-arc camber (fig. 4), 
and a spanwise distribution of twist (fig. 5). In addition to these three basic 
camber elements, two local camber modifications were made. The primary local modifi- 
cation was to add a spanwise bump to reduce the upper-surface curvature; this change 
was added conically. The second local modification was to increase the leading-edge 
camber forward of the " leading-edge -device” hinge line shown in figure 1. This addi- 
tional leading-edge camber varied linearly from a value of 0° at the inboard edge of 
the leading-edge device (43-percent span location) to a value of 5° (positive lead- 
ing edge down) at the wing tip. These camber elements constitute the basic cambered 
wing. An alternate leading edge was designed to be identical to the basic wing, 
except that the local leading-edge camber added at the tip was changed from 5° 
to -2°. 


The final step in the design process was to add a balance housing to the com- 
pleted wing geometry. The balance-housing size was minimized to provide the minimum 
flow distortion to the wing flow field. The balance housing was faired smoothly into 
the wing, both longitudinally and laterally. The final wing design was carried out 
with the balance housing in the computational model. 


WIND-TUNNEL MODEL 

An isolated wing model was sized to fit the Langley Unitary Plan Wind Tunnel. 

The large size of the model helped to achieve surface tolerances of ±0.001 in. on the 
leading edge and ±0.005 in. over the main portion of the model. The wind-tunnel 
model was constructed of aluminum. In table I, the model coordinates for the wing 
with the basic leading edge are given in the format of reference 1 3. The coordinates 
of the model with the alternate leading edge are presented in table II. A steel 
adapter was constructed to affix the internally mounted strain-gage balance to the 
model and orifices for 100 pressure taps were also installed in the model. The loca- 
tions of these pressure taps are listed in table III. 


TEST INFORMATION 

These tests were conducted in the low Mach number test section of the Langley 
Unitary Plan Wind Tunnel, which is a variable Mach number, variable-pressure, 
continuous -flow, supersonic tunnel. The test section is approximately 4.0 ft by 
4.0 ft. (See ref. 10 for a more detailed description of this facility.) Figure 6 
is a photograph of the model installed in the wind tunnel. 
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Tests were conducted at the following nominal test conditions: 


M 

P Q / psf 

P / psf 

oo 

q , psf 

CD 

v ° F 

R, per foot 

1.58 

1072 

260 

454 

125 

2.0 x 10 6 

1.62 

1085 

248 

455 

125 

2.0 x 10 6 

1.66 

1099 

237 

456 

125 

2.0 x 10 6 

1.70 

1113 

226 

456 

125 

2.0 x 10 6 

2.00 

1254 

160 

449 

125 

2.0 x 10 6 


To ensure fully turbulent boundary-layer flow over the model, transition strips 
composed of No. 60 carborundum grit were sprinkled on the upper and lower model sur- 
face 0.4 in. behind the model leading edge (measured streamwise). The transition 
strips were about 0.125 in. wide. The size and location of the transition strips 
were determined by the method of reference 14. 

Angle of attack ranged from approximately 0° to 14°, but most of the pressure 
data were taken between approximately 6° and 14°, inclusive. The measured angle of 
attack was corrected for tunnel-flow angularity and for the deflection of the balance 
and sting under load. Flow-angle corrections were determined by testing the wing in 
both upright and inverted orientations. Pressure data were obtained from six 48-port 
scanning valves mounted outside the tunnel. 

After completing the pressure test, the pressure instrumentation was removed and 
force tests were conducted on the same model. Forces and moments acting on the 
model were measured by means of a six-component strain-gage balance contained within 
the model. The balance was connected through a supporting sting to the model sup- 
port system of the wind tunnel. Two balance -chamber pressure measurements were made 
throughout the force program, and the average of the two chamber pressures was 
applied to the model base area to correct the axial force to a condition of free- 
stream static pressure on the base. After completing the force test, oil-flow photo- 
graphs of the wing upper surface were taken at M = 1.62 for a « 8°, 10°, 12°, and 
14°. 


DISCUSSION OF RESULTS 

The pressure data are discussed first, followed by a discussion of the force and 
moment data. The experimental data used in this discussion are limited to those 
needed for discussion purposes; however, complete plotted and tabulated experimental 
data are presented in appendixes A and B. The associated nonlinear potential-theory 
estimates are for a 57 x 57 grid and a 1-in. marching step. An assessment of the 
effect of grid density and marching-step size on the accuracy and computer execution 
time of the nonlinear potential-theory estimates is the subject of appendix C. 


Pressure Results 


All pressure results are presented as spanwise distributions of pressure coeffi- 
cients. A detailed discussion of the basic leading-edge results is followed by a 
briefer discussion of alternate leading-edge results. 
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Basic Leading Edge 


For the design conditions of a = 12° and M = 1.62, the effects of perturba- 
tions in angle of attack and in Mach number are presented in figures 7 and 8, 
respectively. 

Effect of angle of attack.- Mach 1.62 pressure coefficient results are shown in 
figure 7 for the design angle of attack (»12°) and for angles of attack approximately 
2° below and above the design. Both experimentally measured pressures and theoreti- 
cally predicted (NCOREL) pressures are presented for longitudinal stations of 10.6, 
15.5, 19.9, and 24.4 in. in figures 7(a) to 7(d). Because the theoretically pre- 
dicted pressures represent the goal of the wing-design effort, the quality of the 
agreement between theory and experiment is a validation of the nonlinear potential 
method for this application. 

Both the experimental and theoretical data show that pressures across the entire 
wing are significantly influenced by changes in angle of attack; however, the lower— 
surface pressures exhibit changes only in magnitude, whereas the upper-surface pres- 
sures exhibit changes in both magnitude and in the character of the pressure 
distribution. 

The lower-surface pressure coefficients increase in magnitude with increasing 
angle of attack, as expected, and the quality of the agreement between NCOREL pre- 
dicted values and experimentally measured values is approximately the same for all 
three angles of attack. At the longitudinal station of x = 10.6, the lower— surf ace 
experimental pressure coefficients are somewhat larger than the NCOREL values with a 
maximum error of about 10 percent. However, the agreement at x = 15.5, 19.9, and 
24.4 is virtually identical. At x = 24.4, the most inboard lower-surface pressures 
are predicted higher than the experimental pressures because of a limitation in the 
NCOREL code, which presently must represent the wing wake as a thin, solid-surface 
extension of the trailing edge. 

On the upper surface of the wing, one effect of increasing angle of attack is to 
decrease the pressure, and this effect is most pronounced in the highly nonlinear 
expansion region near the leading edge. Increasing angle of attack can also change 
the character of the pressure distribution, and this is best illustrated by the 
experimental results at x = 19.9 shown in figure 7(c). At the smallest angle of 
attack (a = 9.92°), the pressure distribution shows a well-behaved expansion out- 
board of n » 0.85 followed immediately by a subcritical-type (isentropic) pressure 
recovery inboard. When the angle of attack is increased to a value of 11 .93°, a 
stronger expansion occurs closer to the leading edge, and a constant-pressure plateau 
of supercritical cross flow develops between v\ values of 0.90 and 0.75. On the 
inboard side, the pressure plateau terminates with a rapid pressure recompression; 
this recompression indicates the presence of a cross -flow shock. As the angle of 
attack is further increased to 13.92°, the magnitude of the pressure plateau 
increases, the extent of the plateau increases, and the cross-flow shock moves 
inboard with increased strength. 

The agreement between experimental and predicted (NCOREL) upper-surface pres- 
sures is best in the leading-edge expansion region, with small differences being 
noted for the last two longitudinal stations. At these last two stations, the wing 
leading-edge radii are small, and it is possible that rotational and/or viscous 
effects, which are not accounted for in the nonlinear potential theory, are influenc- 
ing the flow. Additionally, at x = 24.4, the leading-edge expansion peak, which 
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occurs for all three angles of attack, is possibly related to inadequate mesh resolu- 
tion around the leading edge. (See appendix C.) The most notable differences 
between experimental and theoretical upper-surface pressures occur at the cross-flow 
shock, where the potential-flow theory underestimates the cross-flow shock strength. 
This error continues into the subcritical region. During the wing design, it was 
recognized that the isentropic assumptions of the theoretical method would predict 
slightly weaker shock jumps, and this was taken into consideration by imposing more 
stringent limits on the allowable shock strengths. 

At the design angle of attack (a * 12°), the agreement between measured and pre- 
dicted (NCOREL) pressures indicates that the overall design-goal pressure distribu- 
tions were experimentally obtained at all four longitudinal stations. Furthermore, 
this good agreement implies that no flow separation due to either the leading-edge 
expansion or the recompression of the cross flow is present. The oil-flow photo- 
graphs, which are discussed subsequently, also support this view. 

Effect of Mach number .- Experimental pressure coefficient results for four Mach 
numbers at the design angle of attack (ct » 12°) are shown in figure 8 along with 
theoretical (NCOREL) estimates. The experimental data show that the basic nature of 
the flow does not change for perturbations about the design Mach number, and that the 
effects of Mach number are generally confined to the supercritical cross-flow region 
near the leading edge on the upper surface. The magnitude of the expansion pressures 
decrease with increasing Mach number, which is the proper trend. Also, the trends of 
the experimental data are accurately predicted by the theoretical (NCOREL) estimates. 

Linear-theory analysis . - Experimental pressure coefficient data at the design 
condition (M = 1.62 and a = 12°), and at angles of attack 2° above and below the 
design, are repeated in figure 9 along with theoretical pressure-coefficient esti- 
mates from a modified Woodward supersonic linear-theory analysis method (ref. 11) 
which includes thickness effects. Near the leading edge, the large pressure gradi- 
ents and extremely low pressures estimated by linear theory show the dramatic effect 
of the subsonic leading-edge singularity. Also, the linear-theory method cannot be 
used to calculate shocks, so the supercritical-subcritical nature of the upper- 
surface flow is not shown. Comparisons of the experimental and linear-theory pres- 
sures illustrate the inability of linear theory to produce any meaningful information 
on the upper-surface pressure distributions resulting from supercritical cross flow 
about wings. 


Alternate Leading Edge 

The alternate leading edge has less leading-edge camber than the basic leading 
edge, and, as discussed in the section entitled "Aerodynamic Design," the camber 
differences are largest at the wing-tip leading edge. These leading-edge camber 
differences are reflected in the spanwise pressure distributions shown in figure 10. 
In this figure, experimental and theoretical pressures are shown for both leading- 
edge geometries at the design conditions of M = 1.62 and a « 12°. The geometry is 
identical for each leading edge between the wing apex and the x = 10.6 position, 
and this is reflected in the identical pressure distributions of figure 10(a). The 
reduced camber of the alternate leading-edge results in the lower leading-edge expan- 
sion pressures as shown in figures 10(b) to 10(d). In general, the quality of the 
agreement between experiment and theory is the same for the alternate leading edge as 
was previously found for the basic leading edge; the most noticeable difference in 
agreement between experiment and theory occurs at the x = 24.4 station, where the 
large-expansion pressure peak predicted is not experimentally measured. 
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Force and Moment Results 


Basic Leading Edge 

Longitudinal force and moment data are presented in figure 11 for the design 
Mach number of 1.62* In addition to the experimental data, predicted results from 
the nonlinear potential method (NCOREL, ref. 6) and from the linear potential-flow 
method (ref. 11) are also shown. The NCOREL estimates of lift and drag include an 
axial-force contribution due to skin friction (C^ = 0.0069 at M = 1.62), which was 
obtained from the method of reference 15. The skin-friction contribution is assumed 
to be invariant with angle of attack. The linear potential-flow drag estimate is the 
sum of the drag due to lift from the method of reference 11, the far-field wave drag 
obtained for an uncambered wing with the same thickness using the method of refer- 
ence 16, and the skin-friction drag from the method of reference 15. 

The experimental lift and moment data in figure 11 are linear with angle of 
attack through about 9° or 10°. Above this angle of attack, the lift-curve slope and 
the moment-curve slope decrease. In general, the experimental force and moment data 
and the NCOREL calculations agree well; however, small differences between these 
results occur at the higher angles of attack. These differences seem to be traceable 
to the disparity between the calculated and experimentally measured cross-flow shock 
strength; specifically, this disparity would cause an overestimation of the lift and 
a consequent overestimation of the drag and a more nose-down pitching moment, since 
the affected portion of the wing is generally aft of the moment reference point. 

These trends can be seen in figure 1 1 . 

The linear potential-theory estimates are also included in figure 11. The 
linear theory overpredicts C L , C D , and longitudinal stability. The linear-theory 
estimates would be somewhat worse had not the vacuum limit been artificially imposed 
in the computer code. It is informative to relate these linear-theory force and 
moment estimates to the pressure estimates shewn in figure 9; the force and moment 
results are much more accurate than the pressure data might suggest. Also, calcula- 
tion by the nonlinear potential method yields a lower C Q than the linear potential 
method, and the more optimistic nonlinear drag value is supported by the e^qcerimental 
data. 


Figure 11(c) presents the drag polar for the experimental data and the two 
potential theories along with an experimentally -derived polar for the equivalent flat 
plate. The equivalent-flat-plate polar, which is calculated from the equation 
C D = C D + tan (a - (Xq), does not include leading-edge thrust and is taken as the 
lower bound on wing performance. At the design of 0.4, the cambered wing shews a 
21 -percent decrease in drag due to lift compared with the equivalent flat wing. Fig- 
ure 11(c) illustrates that the application of this technology to advanced aircraft 
could provide significant benefits for supersonic maneuvering. Also, the linear- 
theory drag polar is optimistic in the high-lift-coefficient range. 

Since the wing leading edge was rounded, which is in contrast to the sharp lead- 
ing edges of typical supersonic wings, it was suspected that the small-disturbance 
assumptions of the far-field wave-drag prediction method might be violated locally 
and that the calculated wave-drag values should be used with caution. To gain 
further insight into this matter, volumetric wave-drag estimates for an equivalent 
uncambered wing were calculated using the nonlinear potential code (NCOREL), the 
linear-theory near-field method (ref. 11), and the far-field wave-drag method. A 
comparison of the three different wave-drag estimates is shown in figure 12? however, 
since an uncambered version of the cambered wing was not constructed, no experimental 
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data are available. At the design Mach number of 1.62, the far-field wave drag is 
about 20 percent higher than that predicted by the NCOREL code, and this difference 
is reflected in the predicted zero-lift drag values shown in figure 12. The near- 
field wave-drag estimate is totally erroneous, apparently because of an inaccuracy in 
the computation of the longitudinal perturbational velocity component at the leading 
edge of the wing. The NCOREL wave-drag estimates are not affected by Mach cone 
limitations. 

The loss in experimentally measured lift and pitching moment, which was pre- 
viously noted at a « 9° or 10° in the discussion of figure 11, coincides with the 
development of trailing-edge separation which was observed in oil-flow patterns. 
Oil-flow photographs for 8°, 10°, 12°, and 14° angle of attack are shown in fig- 
ure 13. The photograph at a = 8° indicates that smooth, attached flow exists 
everywhere on the wing with the exception of a very small region of separated flow 
at the wing- tip trailing-edge location. At a = 10°, the flow pattern changed 
only slightly, but the separated region on the outboard portion of the wing trailing 
edge enlarged, and a new region of incipient wing trailing-edge separation formed 
inboard. At a = 12°, the smoothly turning flow behind the wing leading edge was 
replaced by a "scalloped" pattern, which possibly indicates the presence of a cross- 
flow shock. At this larger angle of attack, the trailing-edge separation regions 
were enlarged. At a = 1 4° , the scalloped leading-edge pattern moved forward toward 
the wing apex, and virtually the entire trailing edge of the wing separated. 

The onset of trailing-edge flow separation has been correlated with a criterion 
presented in reference 12. This criterion relates the minimum pressure coefficient 
allowable for attached flow at the trailing edge to the free-stream Mach number and 
trailing-edge sweep angle. This trailing-edge criterion is shown in figure 14. The 
experimentally measured plateau pressure coefficient for three angles of attack is 
shewn on the left-hand side of the figure. It is also shown in figure 14 that the 
onset of trailing-edge separation as shown in the oil— flow photographs of figure 13 
correlates well with the empirically determined criterion for the present condition 
of M = 1.62 and a trailing-edge sweep angle of 33°. 


Alternate Leading Edge 

Longitudinal force and moment data are presented in figure 15 for the basic and 
alternate leading-edge configurations at the design Mach number of 1.62. At 12° 
design angle of attack, there is no significant difference in the forces and moments 
produced by the two configurations; a close examination of the tabulated data indi- 
cates that the basic leading-edge configuration has perhaps two counts less drag than 
the alternate leading-edge configuration. 

The most significant difference between the two configurations is shown in fig- 
ure 15(b), where the alternate leading-edge wing produces the lower drag at low lift 
coefficients and produces the higher maximum lift-drag ratio. Both these differences 
are a direct result of the reduced camber drag for the alternate leading-edge config- 
uration compared with the basic leading-edge configuration. 


CONCLUDING REMARKS 

The experimental results of this report represent a verification of a design 
procedure for efficient, high-lift wings at a supersonic design point where Mach 
number is 1.62, angle of attack is 12°, and lift coefficient is 0.4. Efficient high 
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lift is achieved by maintaining attached supercritical cross flow over a major por- 
tion of the wing and then recompressing to subcritical cross-flow conditions through 
a controlled cross-flow shock. This process does not create boundary -layer separa- 
tion. The actual design process, which relies upon nonlinear potential-flow methods, 
is described in detail, and the cortparisons with experimental surface-pressure data 
and longitudinal force and moment data confirm the accuracy of the design method. 

Results are presented which show that design conditions of Mach number and angle 
of attack could be varied slightly without changing the desired flow structure and 
that the nonlinear potential method could accurately predict the change in pressure 
and forces caused by these variations. A disparity between the experimental cross- 
flow shock strength and the calculated isentropic cross-flow shock strength at 
Mach 1.62 is shown in the pressure comparisons, and that disparity produced a small 
overestimation of lift and drag at the higher angles of attack and higher levels of 
longitudinal stability than those measured. Further comparisons of the experimental 
data at Mach 1.62 were made with linear-theory estimated, results. The poor quality 
of the linear potential-theory pressure estimates was noted, but the integrated force 
conparisons were more accurate than the pressure results might indicate. These 
comparisons showed that linear theory is useful as a preliminary performance analysis 
tool but that stability and design studies require a more sophisticated approach for 
the conditions of this study. Oil-flow photographs at Mach 1.62 showed a region of 
trailing-edge separation at high angles of attack, and the experimental pressure data 
were correlated with a trailing-edge separation criterion. This correlation showed 
that the onset of trailing-edge separation was predictable and could be controlled 
through planform, camber surface, angle of attack, Mach number, or a combination of 
these parameters. The overall efficiency of the wing was quantified at the design 
Mach number (1.62) by comparing the experimentally measured drag polar with the 
equivalent flat-plate drag polar (0 percent leading-edge thrust). At the design lift 
coefficient (0.4), the attached-flow, cambered-wing concept showed a 21 -percent per- 
formance improvement relative to the equivalent flat wing. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 25, 1984 
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Figure 11.- Experimental and theoretical longitudinal forces 
and moments for basic leading-edge wing at M = 1.62. 
















Figure 12.- Comparison of zero-lift wave-drag estimation methods. 
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(a) a = 8°. 

Figure 13.- Oil-flow photograph of basic leading-edge wing at 


L-84-50 
M = 1.62. 
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Aft portion 
(b) a = 10°. 
Figure 13.- Continued. 
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(b) a = 10°. 
Figure 13.- Continued. 
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Aft portion 
(c) a = 12°. 
Figure 13.- Continued. 
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(c) a = 12°. 
Figure 13.- Continued. 
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Aft portion 
(d) a = 14°. 
Figure 13.- Continued. 
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(d) a = 14°. 
Figure 13.- Concluded. 


L-84-57 
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Figure 14.- Critical trailing-edge pressure estimates from reference 12. 


56 






(a) C m and a versus C L . 

Figure 15.- Experimental longitudinal forces and moments for 
basic and alternate leading-edge wings at M » 1.62. 










APPENDIX A 


EXPERIMENTAL DATA PLOTS 


The pressure-coefficient data are plotted against the nondimensionalized span- 

V. / T?rn7\ A « rj \ “ 


wise coordinate h (ETA in figures), 
plotted in summary form in figures A1 


Hie entire set of pressure-coefficient data is 
. - and A2. Crossplots of the pressure coefficient 

are shown m figure A3. A summary of the longitudinal force and moment data are 
plotted m figures A4 and A5. 
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(a) Continued. 


Figure A1.- Continued. 
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(a) Continued. 


Figure A1.- Continued. 
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POINT ALPHA 



ETA 

X=24.4 IN 
(a) Concluded. 


Figure At.- Continued 
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Figure A1.- Continued. 
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( d ) Continued . 


Figure A1 Continued. 
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POINT ALPHA 



0 .2 A .6 .8 1.0 

ETA 


X-15.5 IN 

(b) Continued. 
Figure A2.~ Continued. 
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Figure A2*- Continued. 
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Figure A2.- Continued 
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APPENDIX B 


EXPERIMENTAL DATA TABULATION 

The experimental RUN, POINT, and Mach numbers, and the angle-of-attack condi- 
tions are tabulated in table B1* Table B2 contains the pressure-coefficient data 
from the experimental program. The data are listed by POINT number, which indicates 
a unique Mach number and angle of attack for a configuration, that is, basic or 
alternate leading edge. The POINT numbers also appear on the Cp plots in appen- 
dix A and in the main text so that the reader can cross-reference the tabulated data 
with the plotted data. Table B3 contains the longitudinal force and moment data 
listed by RUN number, which indicates a variation of angle of attack at a constant 
Mach number for a configuration. The RUN numbers appear on the force and moment 
plots in appendix A and in the main text so that the reader can cross-reference plot- 
ted and tabulated results. 

POINT 1124 is the corrected POINT 124. Analysis of the experimental data 
revealed that the wind-tunnel operating conditions unexpectedly surged by 3 percent 
while this data point was recorded. The wind-tunnel operating conditions did not 
vary by more than 0.1 percent for the other points in this run; therefore, POINT 124 
was corrected to the average operating condition of the other points in this run. 

Both the original and corrected data for POINT 124 are tabulated, but only POINT 1124 
is plotted. 
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TABLE B1 EXPERIMENTAL RUN SCHEdQTe P° 0R QUALITY 


Basic leading edge 


Run Point Mach Alpha 
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Run Point Mach Alpha 


5 . 

93 

7 . 

92 

8 . 

91 

9 . 

91 

10 . 

91 

11 . 

90 

12 . 

91 

13 . 

91 

5 . 

91 

7 . 

90 

8 . 

93 

5 . 

97 

7 . 

96 

8 . 

97 

9 . 

93 

10 . 

95 

11 . 

93 

12 . 

95 

13 . 

95 

5 . 

98 

5 . 

98 

5 . 

99 

7 . 

99 

9 . 

01 

9 . 

99 

10 . 

97 

11 . 

98 

12 . 

98 

13 . 

98 

5 . 

99 

5 . 

92 

7 . 

94 

8 . 

97 

9 . 

96 

10 . 

94 

11 . 

94 

12 . 

,91 

13 . 

,91 

5 . 

,91 

5 . 

,80 

7 . 

,81 

9 . 

,82 

11 . 

,80 

13 . 

,81 

5 . 

,80 


Reference values adjusted, see preceding page 
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APPENDIX C 


EFFECT OF GRID DENSITY AND STEP SIZE ON NONLINEAR 
POTENTIAL THEORY (NCOREL) RESULTS 

An important consideration for any finite-difference computer program is the 
number of grid points necessary to accurately resolve the given problems at minimum 
cost. In this appendix, the results of a systematic variation of grid-spacing param- 
eters are presented. These results include plots of spanwise pressure distributions 
which compare NCOREL calculations with experimentally obtained data and a table of 
the integrated force and moment coefficients with computer execution times. The 
experimental force and moment data in the table are interpolated to a = 11.9 3 °, and 
the skin-friction axial force of 0.0069 has been removed. 

The grid-spacing parameters assessed in this appendix are the grid density, 
which is held fixed for each two-dimensional cutting plane, and the spherical 
marching-step size DR, which is the distance between each two-dimensional cutting 
plane. The grid density is specified as M x n where M is the number of grid points 
on the body and N is the number of grid points from the body to the outer boundary 
(bow shock). The computational plane grid consists of evenly spaced grid points, but 
in the physical plane the grid points are concentrated near the leading edge to more 
accurately resolve the large leading-edge flow gradients. The NCOREL code marches 
implicitly along spherical cutting planes which are specified at increasing radii 
from the apex of the geometry. The implicit marching technique theoretically allows 
an infinitely large marching step (i.e., no bounds imposed by the CFL criterion), 
and the use of spherical cutting planes allows the code to be used at somewhat lower 
supersonic Mach numbers than would be the case if a Cartesian system were used. 
Without the bounds of the CFL criterion to limit the marching-step size, as is the 
case for ejqplicit marching techniques for hyperbolic flow, the only restriction on 
the marching step is that it must be sufficiently small to accurately model the 
geometry. 

In figure Cl , computed pressures are compared with experimental data at 
a = 11.93° and M = 1.62 for three different grid densities and a 1-in. marching 
step. The increase in grid density from 15 x i 5 to 29 x 29 strongly affects the 
calculated pressure distribution, especially around the leading edge where the gradi- 
ents are strongest. Also, the resolution of the cross-flow shock is quite poor for 
the 15 x is grid. The effect of the increase in grid density from 29 x 29 to 
57 x 57 is not as noticeable on the first two spanwise sections, which are relatively 
thick, but is apparent on the upper surface of the last two sections. In general, 
the effect of increasing the grid density is to provide more accurate spanwise pres- 
sure calculations and a sharper definition of the supercritical and subcritical 
cross-flow regions. 

The effect of three marching- step sizes DR on the NCOREL pressure estimates is 
shown in figure C2 for a constant 57 x 57 grid density. The primary effect is a 
slightly improved cross-flcw shock definition for decreasing step size. At the most 
aft spanwise section, the smallest step size provides the most accurate definition of 
the trailing edge, and the effect on the pressure distribution near the trailing edge 
is apparent. 
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The integrated force and moment estimates are cataloged in table Cl . The most 
ixpensive NCOREL case (DR = 0,5 in., 57 x 57 grid) does not necessarily agree best 
\fith the experimental data. As pointed out in the main body of this paper, this 
jrror is in large measure due to the disparity between the calculated isentropic 
sross-flow shock strength and the experimentally measured cross-flow shock strength, 
^t is important to note that accurate force and moment estimates can be obtained for 
Relatively small run times (DR = 1.0 in., 29 x 29 grid). 
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TABLE Cl SUMMARY OF GRID DENSITY AND STEP-SIZE EFFECTS ON NCOREL 
FORCE AND MOMENT ESTIMATES AND EXECUTION TIMES 

[Experiment interpolated to a = 11.93°; Cf = 0.0069 removed] 


Step size, 
in. 

Grid 

size 

C N 

Ca 

c L 

C D 

c m 

CDC Cyber 175 
CPU seconds 

1 .0 

15 x 

15 

0.40334 

-0.02881 

0.40059 

0.05520 

-0.04148 

134.7 

1 .0 

29 x 

29 

.40300 

-.02751 

.39999 

.05640 

-.03928 

339.5 

1 .0 

57 x 

57 

.40558 

-.02693 

.40238 

.05748 

-.03974 

1881 .8 

.5 

57 x 

57 

.40394 

-.02677 

.40075 

.05731 

-.03887 

3537.0 

2.0 

57 x 

57 

.40660 

-.02711 

.40342 

.05752 

-.04103 

1718.1 

Experiment 



.3929 

-.0285 

.3903 

.0533 

-.0328 
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